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We use the density-functional theory to calculate the total energy of mixed crystals (Ga,Mn)As 
with a small concentration of various donors. We find that the formation energy of Mn depends 
strongly on the partial concentrations of Mn in the substitutional and interstitial positions, and on 
the concentration of other dopants. The composition dependence of the formation energies repre- 
sents an effective feedback mechanism, resulting in the self-compensation property of (Ga,Mn)As. 
We show that the partial concentrations of both substitutional and interstitial Mn increase propor- 
tionally to the total concentration of Mn. 
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I. INTRODUCTION 

The dilute magnetic semiconductors (DMS), such as 
GaAs doped with a large amount of Mn, represent an 
important class of mixed crystals with promising applica- 
tions in spin electronics Q. The ferromagnetic behavior 
of these materials is mediated by the holes in the valence 
band |1, Q] ■ It is sensitive to the number of free carriers 
and to the level of charge compensation. Mn atoms sub- 
stituted in the cation sublattice of a III-V semiconductor 
are acceptors and produce one hole each. It is known, 
however, that some Mn atoms occupy the interstitial po- 
sitions and act as double donors [1, IE 0) S H ■ 

The interplay between the substitutional and intersti- 
tial incorporation of Mn into the GaAs lattice has, to- 
gether with co-doping (cf. Ref. 9]), a crucial effect on 
the physical properties of the mixed crystal (Ga,Mn)As. 
A systematic study of formation energies of the substi- 
tutional (MuQa) and interstitial (Mnjnt) manganese can 
help to understand it on a microscopic level. 

In the case of weak doping small changes in the im- 
purity concentration can easily move the Fermi energy 
Ep across the band gap with a negligible influence on 
the density of states. That is why the dependence of the 
formation energies on the number of electrically active 
impurities is usually represented by their dependence on 
Ep. The Fermi-level dependent formation energy is ob- 
tained by adding (or subtracting) AEp to the formation 
energy calculated for a particular electronic configuration 
|10| . In the case of a strongly doped and mixed crystals, 
however, the redistribution of the electron states in the 
valence band due to the impurities cannot be neglected 
and the density-of-states effect modifies the simple Fermi- 
level rule for the formation energies. 

That is why the knowledge of the formation energy 
as a function of the impurity concentrations is neces- 
sary in the DMS. To calculate the energy needed to 
incorporate Mn and other impurities in a mixed crys- 
tal, we use the trick relating them to the composition- 
dependent total energy of the mixed crystal jTl'|. This 



quantity is obtained within the density-functional the- 
ory for a series of (Ga,Mn)As mixed crystals with vari- 
ous content of Mn in substitutional and interstitial posi- 
tions, and with variable concentration of the compensat- 
ing donors. The use of the coherent potential approxima- 
tion (CPA) combined with the tight-bin ding linearized- 
muffin-tin-orbital method (TB-LMTO) [T^ makes pos- 
sible to change the chemical composition continuously. 
The lattice relaxation around the impurities and the clus- 
tering of the Mn atoms are omitted within the CPA. For 
simplicity, we consider only the interstitial Mn atoms in 
the T(As4) position; the energy of the other T(Ga4) po- 
sition is almost the same [iJl • The formation energies are 
obtained as the first derivatives of the total energy with 
respect to the corresponding partial concentration pd| . 

Assuming a quasi-equilibrium deposition conditions, 
characterized by an effective growth temperature, we use 
the calculated formation energies to estimate the num- 
bers of MuQa and Mni„t in (Ga,Mn) As mixed crystal. We 
also present a simple way to determine the partial con- 
centrations directly from the composition dependence of 
the formation energies, without solving thermodynamical 
balance equations. 



II. COMPOSITION DEPENDENCE OF THE 
FORMATION ENERGIES 

We consider an impure or mixed crystal with several 
kinds of impurities Ii, I2, etc.. The total energy of the 
mixed crystal W{xi, X2, ■■), normalized to a unit cell, de- 
pends on their molar concentration Xi. As we showed 
recently [ij, the formation energy Ei of an impurity I^ 
can be obtained by differentiating W{xi,X2, ■■) with re- 
spect to Xi, namely 



Ei{xi,X2,.-) 



dW{xi,X2, ■ 
dxi 



.^«to™(I.)+^«to™(host), 

(1) 

The last two terms in Eq. (1) are the total energies of an 
atom li and of the corresponding atom of the host, which 
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Figure 1: Relative formation energy AE of the substitutional 
MuGa in Gao.96Mno.04As as a function of the concentration 
of various donors. 
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Figure 2: Relative formation energy AE of the interstitial 
Mn in Gao.96Mno.04As as a function of the concentration of 
various donors. 



has been replaced by I^. 

Generally, the definition of the formation energy is not 
unique and depends on the way in which the atomic ener- 
gies £;"*°'"(host) and £;"*°'"(I,) in Eq. (1) are obtained. 
We use the energies of neutral atoms in their ground 
state. It is important to notice, however, that the addi- 
tional constant in Eq. (1) does not depend on the actual 
chemical composition of the material. It is not important 
for the concentration-dependent trends we have in mind. 

That is why we consider now only the relative forma- 
tion energies AEi, obtained from their actual values Ei 
by subtracting the corresponding formation energy cal- 
culated for the reference material. As a reference, we 
take Gao.96Mno.04As with all Mn atoms in regular MuQa 
positions. 

The composition dependence of the formation energies 
Ei is characterized by the coefficients of the linear expan- 
sion around the reference point, 



K,,. 



dxj 



dxidxi 



(2) 



The correlation energies Kij |ll| form a symmetric ma- 
trix. Kij < means that the presence of the defects I; 
supports formation of Ij and vice versa; Kij > indicates 
the opposite tendency. 

The dependence of the relative formation energy AE 
of the MuGa on the concentration of various donors is 
summarized in Fig. 1. We considered four representative 
examples. Sbas and Sica a-re typical donors with one 
extra electron, situated at anion and cation sublattice, 
respectively. The other two cases, i.e., As antisite defect 
AsGa and Mnint are the most important native defects 
in (Ga,Mn)As, both acting as double donors. Fig. 1 
shows that the formation energy of MuQa decreases in 
the presence of an increasing number of donors. The 
curves are grouped into pairs according to the charge 
state of the donors, with only a minor influence of the 
particular chemical origin of the defect. The dependence 
is almost linear for low concentrations and the slope of 



the function is roughly proportional to the charge state 
of the donor. This all indicates that the variations of 
the formation energy of Mnca are mostly determined by 
the Fermi-level effect, not by the redistribution of the 
density of states induced by the other defects. 

Analogous results are obtained for the formation en- 
ergy of the interstitial Mn in the T(As4) position, as 
shown in Fig. 2. In this case, however, the formation 
energy of Mnjnt increases with increasing number of the 
donors. It is important to notice that the steep increase 
of the formation energy represents a feedback mechanism 
limiting efficiently the number of Muint- The same is 
valid also for the formation energy of the As antisite de- 
fect. 

Fig. 3 shows the formation energies of the two native 
defects, i.e. Asca and Mnint , in the mixed crystal with 
varying number of the substitutional Mn in the Ga sub- 
lattice. Notice that both relative quantities, being pinned 
to zero for the reference material with 4 % of MuQa, are 
almost identical. In both cases, the formation energy is a 
decreasing function, indicating an increasing probability 
of formation of these defects in materials with a higher 
concentration of Mn. 

This self-compensation tendency is a very important 
mechanism controlling the basic physical properties of 
(Ga,Mn)As mixed crystals. It is the reason for the ob- 
served low doping efficiency of Mn in GaAs [lj|. The 
increasing number of both Asca and Mujnt also explains 
the expansion of the lattice of (Ga,Mn)As with an in- 
creasing concentration of Mn M^ . 



III. DYNAMICAL EQUILIBRIUM BETWEEN 

Mnca AND Mnint 

In this Section, we use the calculated formation en- 
ergies to simulate the incorporation of Mn into the 
(Ga,Mn)As mixed crystal. We assume that the probabil- 
ities that an Mn atom occupies either substitutional or 
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Figure 3: The formation energy of the interstitial Mn and As 
antisite defect in Gai-ajMn^iAs as a function of the concen- 
tration X of Mn atoms substituted in the Ga sublattice. The 
formation energies are referred to the values corresponding to 
Gao.96Mno.04As 
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Figure 4: The partial concentrations of Mnca (solid) and 
Muint (dashed) as a function of the total concentration of 
Mn in (Ga,Mn)As for T^// = 500 K. The influence of other 
dopants is not considered. 



interstitial position are related to the corresponding for- 
mation energies Es and Ej also in the non-equilibrium 
epitaxial growth. As a simplest approximation, we char- 
acterize the deposition condition with some effective 
temperature T^ff and use the corresponding Boltzmann 
weighting factors. 

To this purpose, the absolute formation energies of 
Mn in the two crystallographic positions, i.e. MuQa a-nd 
Muint, are required. They are obtained for the reference 
system Gao.96Mno.04As including the additive terms from 
Eq. (1). For the epitaxial growth, it is reasonable to use 
£jatoTn calculated for isolated atoms as stated above. A 
linear interpolation for the dependence of Es and Ej on 
the corresponding partial concentrations xs and xj is 
used: 



Es{xs,xi) 



E°s 



+ ussxs + Ksixi 



Ei{xs,xi) = E°^ + Kisxs + Kiixi 



(3) 



(4) 



with E"^ 



0.31 eV, E° = 0.42 eV, K55 = -0.17 eV, 
Ksi = f^si = —6.03 eV, and kjj = 10.33 eV. 

The composition dependent formation energies, Eqs. 
(3-4), define the thermodynamic probabilities ps and pj 
that extra Mn atoms occupy substitutional or interstitial 
position in the mixed crystal with a given composition. 
They are 



PS, I 



e-xp{-Esj/kTeff) 



eM-Es/kT,ff) + eM-Ei/kTeff) ' 



(5) 



On the other hand, these probabilities determine the 
number of Mn atoms that substitute for Ga or occupy 
the interstitial positions. The resulting changes of the 
partial concentrations xs and xi due to the variation dx 
of the total concentration x of Mn are 



and the dependence of xs and xi on x can be obtained 
by the integration of Eq. (6). 

Fig. 4 shows the solution of Eq.(6) for Te// = 
500 K. For the lowest concentrations {x < 0.015), Mn 
atoms occupy preferentially the substitutional positions 
which have lower formation energy (cf. Eqs. (3,4)). For 
higher concentration of Mn, however, the difference of 
Ej{xs,xj) — Es{xs,xi) decreases and approaches zero. 
From this point, both positions can be occupied with a 
comparable probability, and the partial concentrations 
of both MuGa and Mnjnt increase proportionally to x. 
Any deviation from the situation with Es = Ej changes 
the formation energies of Mnca and Mnjnt in such a way 
that the dynamical equilibrium is restored. As a result, 
the high-concentration regime with co-existing MuQa and 
Muint is stabilized. This finding does not depend much 
on Teff over a wide temperature range. We can conclude 
that the partial concentrations of Mn are simply given by 
the equation 



Esixs,xi) = Ei{xs,xi), 



(7) 



dxf 



Psdx, dxi = pidx, 



(6) 



together with the condition xs + xj = x. Combining 
Eqs. (3,4) with (7) we find that the proportion of MuQa 
and Mujnt atoms is roughly 3:1, in a very good agreement 
with both experiment ^ and theoretical expectations [6|, 
■15j for the as-grown materials. 

It is important to point out that according to Eqs. 
(3,4,7) both MuQa and Mujnt atoms remain metastable 
with the activation energy w 0.3 eV in the whole con- 
centration range shown in Fig. 4. The formation ener- 
gies used in our dynamical-equilibrium approach control 
the preferential incorporation of Mn atoms during the 
growth. The annealing process, on the other hand, de- 
pends on the barriers preventing the Mn atoms to leave 
their metastable positions. The barriers are lower for 
Muint than for Mnca position |3|j so that the post- 
growth treatment can substantially reduce the number of 



the interstitial Mn atoms without a remarkable change 
of the number of MuQa- 



IV. SUMMARY 

We have shown that the formation energies of Mn 
in either substitutional or interstitial position depend 
strongly on the partial concentrations of both MuQa and 
Muint, and also on the number of compensating donors. 
Also the formation energy of Asca antisite, the main na- 
tive defect in (Ga,Mn)As, is very sensitive to the concen- 
tration of Mn. 

The composition dependence of the formation ener- 
gies represent a feedback mechanism which defines a dy- 
namical equilibrium between Mnca, Mnint, and other de- 



fects and impurities during the growth. In particular, 
we found that at higher Mn concentrations the number 
of both Mnca and Muint increases proportionally to the 
total concentration of Mn in the as-grown (Ga,Mn)As 
mixed crystal. 

In addition, the concentration dependence of the for- 
mation energy of the AsQa antisite defects indicates that 
an increasing number of these donors also participate 
in the compensation of the regular Mnca acceptors for 
higher Mn concentrations. 

Acknowledgment This work has been done within the 
project AVOZl-010-914 of the AS CR. The financial sup- 
port was provided by the Academy of Sciences of the 
Czech Repubhc (Grant No. A1010214) and by the Grant 
Agency of the Czech Repubhc (202/04/583). 



[1] H. Ohno, J. Crystal Growth 251, 285 (2003). [10] 

[2] T. Dietl, A. Haury, Y. Merle, Phys. Rev. B 55, R 3347 

(1997). [11] 

[3] T. Dietl, H. Ohno, F. Matsukura, J. Cibert, D. Ferrand, 

Science 287, 1019 (2000). [12] 

[4] K. M. Yu, W. Walukiewicz, T. Wojtowicz, I. Kuryliszyn, 

X. Liu, Y. Sasaki, J. K. Furdyna, Phys. Rev. B 65, R 

201303 (2002). [13] 

[5] J. Masek, F. Maca, Acta Phys. Pol. A 100, 319 (2001). [14] 

[6] F. Maca, J. Masek, Phys. Rev. B 65, 235209 (2002). 
[7] S. Sanvito, G. Theurich, N. Hill, J. Superconductivity 

15, 85 (2002). [15] 

[8] J. Blinowski, P. Kacman, Phys. Rev. B 67, 121204 

(2003). [16] 

[9] K .M. Yu, W. Walukiewicz, T. Wojtowicz, W .L. Lim, 

X. Liu, U. Bindley, M. Dobrowolska, J .K. Furdyna, Phys. 

Rev. B 68, R41308 (2003). 



P. Mahadevan, A. Zunger, Phys. Rev. B 68, 075202 
(2003). 

J. Masek, I. Turek, V. Drchal, J. Kudrnovsky, F. Maca, 
Acta Phys. Pol. A 102, 673 (2002). 

I. Turek, V. Drchal, J. Kudrnovsky, M. Sob, P. Wein- 
berger, Electronic Structure of Disordered Alloys, Sur- 
faces and Interfaces, (Kluwer, Boston 1997). 
J. Masek, F. Maca, Phys. Rev. B 69, 165212 (2004). 
B. Beschoten, P. A. Crowell, I. Malajovich, 
D. D. Awschalom, F. Matsukura, A. Shen, H. Ohno, 
Phys. Rev. Lett. 83, 3073 (1999). 

J. Masek, J. Kudrnovsky, F. Maca, Phys. Rev. B 67, 
153203 (2003). 

S. C. Erwin, A. G. Petukhov, Phys. Rev. Lett. 89, 227201 
(2002). 



